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Рассмотрен общий принцип инвариантности в теории переноса излучения. Получены инвариантные соотношения, справедли- 
вые для дисперсных пространственно-ограниченных сред, которые обобщены на случай среды, ограниченной отражающими 
поверхностями. Определен коэффициент диффузного отражения излучения от дисперсной среды с использованием принципа 
инвариантности. 


Принципы инвариантности и симметрии игра- 
ют важную роль в современной физике [1]. В опти- 
ке инвариантные соотношения широко использу- 
ются при расчете оптических систем [2]. В теории 
переноса излучения принцип инвариантности 
применяется для получения результатов даже без 
использования уравнения переноса излучения 
[3, 4]. Общий принцип инвариантности формули- 
руется следующим образом [5]. Пусть имеется рас- 
сеивающая среда с произвольно распределенными 
источниками, мощность которых зависит только 
от одной пространственной координаты г (т - оп- 
тическая плотность среды). Выделим внутри среды 
два уровня ту и т 2 и рассмотрим слой, лежащий 
между ними. Тогда интенсивность диффузного из- 
лучения (или функция источника) в этом слое од- 
новременно есть функция источника всей среды. 
Излучение на границах выделенной области опре- 
деляется рассеянием в остальной части дисперсной 
среды. Частными случаями этого принципа явля- 
ются классические принципы инвариантности 
В.А. Амбарцумяна [6] и С. Чандрасекара [7] или 
метод удвоения слоев, используемый в методе мно- 
гократных отражений [8]. Важным моментом об- 
щего принципа инвариантности является то, что 
он применим для среды любой геометрии, т.к. вы- 
деленная область, являющаяся частью всей среды, 
может быть ограничена поверхностью произволь- 
ной формы, и эта часть взаимодействует со всей 
средой благодаря поступлению излучения через 
границы выделенной части. Эта же идея содержит- 
ся в методе поверхностных псевдоисточников [9]. 
Однако это свойство общего принципа инвариант- 


ности не использовалось для получения инвариа- 
нтных соотношений в ограниченных средах. 

Известные инвариантные соотношения форму- 
лируются для среды, неограниченной в попереч- 
ном (по отношению к направлению распростране- 
ния излучения) направлении. В данной работе по- 
лучено соотношение, которое остается постоян- 
ным при изменении поперечных оптических раз- 
меров среды и коэффициента отражения поверх- 
ностей, ограничивающих рассеивающую среду. 
Наиболее близким к найденному инварианту явля- 
ется фотометрический инвариант [10, 11], имею- 
щий вид 


1 + р 2 +г 2 
— — г = сопві:, 


( 1 ) 


где г и р - коэффициенты пропускания и отраже- 
ния слоя дисперсной среды. Рассмотрим, как вли- 
яет поперечная ограниченность среды, располо- 
женной между отражающими поверхностями с ко- 
эффициентом отражения г, на пропускание и отра- 
жение излучения. 

Введем основные параметры, поясняющие пос- 
тановку задачи. Используем модель среды в виде 
прямоугольного параллелепипеда с оптическими 
размерами т ѵ т у> , г,., освещаемого коллимирован- 
ным потоком излучения интенсивности / 0 = 1 , нап- 
равленным по нормали к плоскости у%, при ис- 
пользовании декартовой системы координат. В ка- 
честве параметров элементарного объема среды 
рассматриваются коэффициент ослабления 
а=(7+ ае (где сх и ае - коэффициенты рассеяния и 


41 


Известия Томского политехнического университета. 2003. Т. 306. № 3 


поглощения), вероятность выживания кванта 
Л= а! а и индикатриса рассеяния излучения х(Ѳ) 
( Ѳ - угол рассеяния), представленная в виде интег- 
ральных параметров 77, Д /л с условием нормировки 
г/ +Р + 4/и=1. Параметры 77, Д /и определяются как: 


Ч = 


к 
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К 
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, где 


5, (т у ,Ѳ,г) = I - 
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К' 1 (т у ,Ѳ,г)т у (А 1 - Д)ехр[-Д (т у ,Ѳ,г)т ѣ ]} 

(1-г)[1 + Д(Ѳ)]{1+ехр[-<(Ѳ)г Ѵ|| ]} 
Ѵ ’ Г {Д+Де Х р[-^(Ѳ)г Ѵо ]}<(Ѳ)г Ѵо 
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Формулы для расчета потоков излучения, пропу- 
щенного /|(го,А,Ѳ/) и отраженного / 2 (то,Л, Ѳ,г) средой 
с отражающими поверхностями, получены с исполь- 
зованием принципа инвариантности и имеют вид 

І х {т 0 ,Р,Ѳ,г) = 

_ ^-^ 2 (т У|г ,А,б,^)еxр[-^(г ; ,_,Л,е,^)г x ^ 

П 2 -5 2 ехр[-2іГ( г ѵ 2 , Л, Ѳ, г )г Іо ] 

і 2 (т 0 ,А,Ѳ,г) = 

і 2 (т 0 ,А,Ѳ,г) = 

і[^у, г Л,Ѳ,г)А-В]ехр[-К(т у ^А,Ѳ,г)т Хо ] 
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где А =\-гК(т к ,А,Ѳ,г), В = К(т ул ,А,Ѳ,г)-г. 

Здесь г 0 - оптические размеры среды, (ж г- коэф- 
фициенты пропускания и отражения отражающих 
поверхностей, К(т Уі „А,Ѳ,г) - функция, определяю- 
щая перераспределение излучения с направлений у, г, 
на направление х, К(т„,А,Ѳ,г) - функция, имеющая 
смысл коэффициента отражения среды при г->оо. 

В частном случае консервативной (Л=1) сим- 
метричной среды (г,:=г..), ограниченной отражаю- 
щими поверхностями с одинаковым коэффициен- 
том отражения г, параметры, входящие в формулы 
(2) и (3), имеют следующий вид: 


ѵ(т у ,Ѳ,г) = 


8Д (т у ,Ѳ,г)8 2 (т у ,Ѳ,г) 
[\-5 х (т у ,Ѳ,г)][2-5 2 (г у ,Ѳ,г)] + ' 


+8 8 х (т у ,Ѳ,г) 


где 


введены обозначения: А, =1-гЯ,(т у ,Ѳ,г), В, = В,( т у , 
Ѳ,г)-г, А[ =1-/ ‘К[(Ѳ), В[ =\-гК[(Ѳ). 

Переменные коэффициенты определяются по 
формулам 

К\(Ѳ) = 2уііі(\-п + РУ, 

К І (т у ,Ѳ) = ^Р-^4-ѵ 1 (т у ,Ѳ); 

К(г у ,Ѳ) = К[(Ѳ)^-ѵ(т у ,Ѳ)- 


ѵ і (т у ,Ѳ) = 1- 
Я(т у ,Ѳ) = 


([1 + Д(^)]{1-ехр[-<(^)г Уо ]} 

К ( ѳ ) т у„ {1 + К (Ѳ) ехр[-< (Ѳ)т ѣ ] } ' 
\-г] + /3- К(т у ,Ѳ) 


(4) 


І-Г/ + Р + К(т у ,ѲУ 

Параметры В4т у ,Ѳ) и К[(Ѳ) вычисляются по 
формуле (4) при подстановке К { ( т у ,Ѳ) и К\{ Ѳ) соот- 
ветственно. Уравнение для К(т ух ,А,Ѳ,г) и 
К(т ух ,А,Ѳ,г) имеют следующий вид: 

К(г уг ,А,Ѳ,г) = 

= Ѵ4 А 7[1-ѵ(г > . 2 ,Л,Ѳ,г)][1-(7 7 -Д)]; (5) 


К(г уг ,А,Ѳ,г) = 


К(т у г , А, Ѳ, г) - 4/л[1 -ѵ(т у . , А, Ѳ, г)] 




К( г Ѵ 2 ,А,Ѳ,г) + 4ц[1 - ѵ(г у 7 , Л, Ѳ, г)] ' 


( 6 ) 


у г - 


Совместное решение уравнений для К(т у „А,Ѳ,г) 
и К(т ух ,А,Ѳ,г) приводит к соотношению: 

Р + 2 ц = — ѵ 


■(1 -А) .(7) 


2Л 

К (г у 2 ; Л; Ѳ; г)[1 + К(т у 7 ; Л; 0; г) 
1-К(т-А;Ѳ;г) 


В случае консервативной среды Л=1: 

„ _ 1 К(т у/ ,А;Ѳ;г)[1 + К(т у/ ,А;Ѳ;г) 

Р+ И 2 1 -К(т у/ ,А;Ѳ;г) ' (8) 

Так как сумма Р+2 /и постоянна, то эти выраже- 
ния являются инвариантами для данного вещества, 
справедливыми при любых оптических размерах 
рассеивающей среды. 
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Рис. 1. Зависимость параметров й (1) и К (2) от поперечных 
оптических размеров среды 



Рис. 2. Зависимость коэффициента отражения от попереч- 
ного оптического сечения среды: г х : 1) 1СР, 2) 70', 
3) ІО 2 , 4) ІО 4 , а=1, Л=7 



Физический смысл величин, входящих в инва- 
риантные соотношения, заключается в следую- 
щем. Параметр Щт у „А,Ѳ,г) имеет смысл коэффи- 
циента отражения от полубесконечного слоя [12]. 
В случае среды, ограниченной в поперечных нап- 
равлениях, данный параметр - величина не посто- 
янная, а является сложной функцией оптических 
свойств и размеров среды, что иллюстрируется 
рис. 1. К(т ух ,А,Ѳ,г) представляет собой функцию 
ослабления, зависящую от поперечных оптических 
размеров и свойств среды и учитывающую допол- 
нительное приращение интенсивности рассеянно- 
го излучения по сравнению с величиной интенсив- 
ности, определяемой законом Бугера [12]. 

Использование инвариантных соотношений (7) 
и (8) позволяет получить в отличие от широко при- 
меняемых в теории переноса излучения численных 
методов, например, метода Монте-Карло [13], ана- 
литическое выражение для коэффициента отраже- 
ния излучения. Согласно (3), поведение коэффи- 
циента отражения р дисперсной среды при изме- 
нении ее параметров в значительной степени опре- 
деляется функцией К: 


к (* у /А;Ѳ;г) = 

■ _ №,;*;»;г) (а+1н мм ^ 
= 1+ КЩ^ (а + 1)+ (А -А Х-Ц - 


где а = 


'7+ 2ц 
Р+2р' 


При отражении излучения от полубесконечной 
среды К(т ух ,А,Ѳ,г) становится равным коэффици- 
енту отражения всей среды р. Формула (9) подобна 
закону Френеля, описывающему отражение нор- 
мально падающего излучения от границы раздела 
двух сред. 

Вид зависимости коэффициента отражения 
слоя фиксированной толщины т х от поперечного 
оптического сечения среды приведен на рис. 2. 

Таким образом, использование общего принци- 
па инвариантности позволяет получить инвариант- 
ные соотношения (7, 8) для пространственно-огра- 
ниченных сред с отражающими границами. Фото- 
метрический инвариант (1), справедливый при лю- 
бой толщине слоя среды, и не являющийся инва- 
риантом при изменении поперечных оптических 
размеров рассеивающей среды, дополнен инвари- 
антом (7) или (8), сохраняющимся при любых из- 
менениях оптических размеров рассеивающей сре- 
ды. Кроме того, использование этих инвариантных 
соотношений позволило определить коэффициент 
диффузного отражения излучения от дисперсной 
среды. 
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